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Abstract 

We consider a transient random walk on Z d which is asymptotically stable, 
without centering, in a sense which allows different norming for each compo- 
nent. The paper is devoted to the asymptotics of the probability of the first 
return to the origin of such a random walk at time n. 
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1. Introduction 

Let S n , n > 0, be a random walk in Z d generated by independent iden- 
tically distributed steps £ n = (£ n i> • • • > £nd)j n > 1, that is, So = 0, S n = 
£i + . . . + £ n . Denote r = and recursively r n+ i = min{fc > r n : S& = 0}; 
by standard convention min0 = oo. Then r = T\ is the first return to the 
origin of the random walk S n . 

In this paper we study the asymptotic behaviour of 

p n := P{r = n} = P{the first return to zero occurs at time n} 
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as n -> oo. Put G(B) = P{r G 5} and 



V 



= P{r < 00} = G[l,oo) < 1. 



n=l 



The measure w n := P{S" n = 0} on Z + is actually the renewal measure gener- 
ated by the r's, that is, 



u n = Y,nr k = n} = Y,G*\n). (1) 

ifc=0 k=0 

Then 

00 00 00 

J2ns n = 0} = ^G* fc [o,oo) = = _, (2) 

n=0 fc=0 fc=0 ^ 

which implies 

i+E:=i^n = o}- lJ 

The random walk S n is called aperiodic if Z d is a minimal lattice for S* n 
in the sense that, for every e > 0, 

sup |Ee i(A ' 6) | < 1. 

Ae[-7r,7r] d \[-e,£]' 1 

Aperiodicity is clearly no essential restriction, as the state space can always 
be redefined, if necessary, so as to make a random walk aperiodic. 

In the sequel we will be studying aperiodic random walks on Z d which 
are asymptotically stable in the following sense: there is sequence c n = 
(c n i, . . . ,c nd ) such that 

X n := (S n i/c nl , . . . ,S nd /c nd ) 4y = (Y h ...,Y d ), 

where Y is a strictly d- dimensional stable random variable. Since this implies 
that each component of X n is asymptotically stable, we know that each c nr 
is in the class RV(l/a r ) of regularly varying at infinity with index l/a r 
sequences (see, e.g. |3|, Section 1.9]), where a r G (0,2] is the index of the 
univariate stable random variable Y r . Thus C n := Y\i c nr is in RV(r)), where 
7] = Y^i^/ a r > d/2. We need the following local limit theorem, in which g 
denotes the density function of Y. 
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Theorem 1. If S n is an aperiodic random walk on 7L d which is asymptoti- 
cally stable in the above sense, it holds that uniformly for x G Z d 

C n F{S n = x} = g(x 1 /c nl , x d /c nd ) + o(l) as n -)• oo. 

In particular, u n = F{S n = 0} ~ g(0, . . . ,0)/C n as n — )■ oo. 



For d = 1 this is the classical local limit theorem of Gnedenko, see 12 



§ 50]; for the case d = 2 it is proved in [5J, and as remarked there, the proof 
extends in a straightforward way to the case d > 2. It can also be viewed as 
a special case of Theorem 6.4 in jsj], where the more general case of matrix 
norming is treated. 

Since it is known that g and its derivatives are bounded, we deduce the 
following 

Corollary 2. If S n is an aperiodic random walk on Z d which is asymptoti- 
cally stable, there exists a constant k such that 

V{S n = x} < k/C n forallxeZ d , 

and, for every fixed k, 

¥{S n ^ k = x} = F{S n = x} + o(l/C n ) 

as n —¥ oo uniformly for x <E 7* d . 

Remark 3. Since every random walk is transient when d > 3, a result which 
for the simplest symmetric random walk goes back to Polya J 13] , the require- 
ment of transience only features for d = 1 and d = 2, when under our 
assumptions it is equivalent to Y^n=i V^n < 00 ■ 

Suppose we know that P{r = n} G RV(— 7) for some 7 > 1. Then 
G/p is so-called locally subexponential distribution. In this case, the local 
asymptotics of the defective renewal function is described in [l], Proposition 
12]: the following limit exists: 



hm E y, (nj =E^ 



2 ■ 



G{n) (1-p) 

so that P{5* ra = 0} ~ P{r = n}/(l — p) 2 . This is the intuition behind the 
following, which is our main result. 
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Theorem 4. Let S n be an aperiodic, transient random walk on 7L , d > 1, 
which is asymptotically stable in the above sense. Then as n — )■ oo, 



p{ r = n}~(l-p) 2 P{S n = 0} 



2ttd r c m (l-p)M0,..-,0) 



Remark 5. 27ie special case where d > 3,E£ n = and E£ n j£ n j = .By < oo 
lozt/i det i? 7^ reads 

P {T = „}-(!- p) »P {4 . = o}~^^ r -^. 

This was proved in an unpublished communication by one of us, and is quoted 
in Chapter A. 6 offyj; this illustrates the increasing importance of local results 
such as this in Mathematical Physics. To the best of our knowledge, Theorem 
[7] was not proved in the literature even in the case of the simplest symmetric 
random walk on Z d with d > 3. 

In the next section we prove Theorem H] by analytic means via a Banach 
algebra technique; this is the method that was used in proving the above 
special case. Then in Section [3] we give a probabilistic proof capturing the 
most probable way that large values of r occur. 

Note also that the recurrent one-dimensional case d — 1 was first studied 
by Kesten in [lo[ where it was proved in Theorems 7 and 8 that when E£ =0 



and E£ 2 < oo, we have P{r = n} ~ y~^f n 3 ^ 2 as n -)■ oo; the recurrent 
case of convergence to a stable law with index a G [1,2] was also considered. 
Two different approaches for proving this equivalence may be found in 0, 
Theorem 1.2]. 

In dimension 2 Jain and Pruitt 0, Theorem 4.1] proved, assuming zero 
mean and finite covariance B, that P{r = n} ~ 27rv det Bn~ x log - n as 
n — > oo. 

The only local result for dimensions 3 and higher we found is one by 



Kesten and Spitzer [111, Theorem lb] where they proved that P{r = n + 1} ~ 
P{r = n} as n — > oo given E£ = 0. 

In conclusion note that different mechanisms are involved in formation 
of large deviations of r in dimensions d — 1, d — 2 and d > 3. In the 
one dimensional recurrent case, the equation ([1]) says that we deal with the 
renewal process generated by r's where r has an infinite mean. In principle 
the same is true for the case d = 2 but here the tail of r is very heavy, it is 
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slowly varying at infinity. In the case d > 3 transience holds, so that we have 
the renewal process generated by a defective distribution which yields that 
the renewal atoms u n are proportional to p n \ this is the main topic addressed 
to in the present article. 



2. Banach algebra approach 

The proof of Theorem [5] is straightforward, if we assume additionally that 
Y^n=i u n < 1 which is equivalent to p < 1/2. The discrete renewal relation 
dlj) between the w's and p's implies that 

1 

l + u(s) 



l-p(s) 

where we put u(s) = Y^n°=i u nS n and p(s) = Y^=oVn,s n , so that 
This is equivalent to 

n 

Vn = E(- 1 ) fe+1 < (fc) 
k=l 

where is the A;-fold convolution of {«„}. Note that u n is regularly varying 
at infinity so that the defective distribution {u n } is locally subexponential, 
see |l| or p, Chapter 4]. Then, as follows from [6|, Theorem 4.30], 

n fc=i n k=i \ < \ jj 

Note 1 + it(l) = 1/(1 — p), so p n ~ (1 — p) 2- u„ as n — > oo. 

What happens if X^^Li u n > 1? In this case we can't expand u(s)/(l + 
u(s)) as a power series in u(s). Nevertheless this is an analytic function 
in u(s), for all (complex) s in \s\ < 1, because w(l) < oo by the assumed 
transience. We can then apply Theorem 1 from the paper [4j and get the 
same result; this reference is based on Banach algebra techniques. 
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3. Probabilistic approach 

The starting point of the probabilistic proof of Theorem H] is the following 
result, which holds in any dimension. 

Lemma 6. Let a sequence Q n which is regularly varying at infinity be such 
that 

F{S n = x} < Q n for all x G Z d , (4) 

oo 

^Qn < oo, (5) 

n=l 

and, for every fixed k, 

¥{Sn-k = x}= F{S n = x}+ o(Q n ) (6) 

as n — >■ oo uniformly in allx G Z d . Let r n be any fixed unboundedly increasing 
sequence. Then 

F{S n = x, t > n} = (1 - p)F{S n = x} + o{Q n ) 

as n oo uniformly in x G Hj d such that > r n . 

Proof. It is equivalent to prove the relation 

F{S n = x,r<n} = P F{S n = x} + o{Q n ). (7) 

We start with the following decomposition: 

n-l 



F{S n = x,r <n} = ^F{S n = x,t = k} 

fe=i 

n-l 

= ^P{S' n _ fc = x}P{r = A;}. 



fc=i 

For every fixed N, we have 

n— N n-N 



J2nS n -k = x}F{r = k} < ^2F{S n - k = x}F{S k = 0} 

k=N 
n-N 

< QkQn-k, (8) 



k=N k=N 

n-N 



k=N 
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by the condition (]3J). For every fixed k, F{Sk = x} — > as \\x\\ —> oo. 
Together with (J4J) it implies that, for every fixed N, 



n— 1 n— 1 



^ P{5* n _ fc = a;}P{r = k} < ^ F{S n ^ = x}Q k 

k=n-N k=n-N 

n-1 

< Yl Qko(l) = o(Q n ) (9) 

k=n—N 

as n — > oo uniformly in \\x\\ > r n . Finally by for every fixed N, 

N N 

nSn-k = x}F{r = k} = F{S n = x}J2 = k } + °(Qn) (10) 

k=l k=l 

as n — y oo uniformly in all x. 

Combining (THT) — (TTUT) we obtain that 



-. 71—1 

limsup — |^P{S n = x,r = k}-F{S n = x}F{r < N}\ 
' ^ n k=i 

^ n-N 

< limsup — Y QkQn-k- 



7i->oo Qn k v 



By regular variation of Q n at infinity, taking into account that the series (jSJ) 
converges, we conclude that the right hand side can be made as small as we 
please by choosing N sufficiently large; this is also a well-known property 
in the theory of locally subexponential distributions, see e.g @, Chapter 4]. 
Now the proof of (JTj) follows by letting N — > oo. 

Lemma 7. Under the conditions of Lemma® as n — > oo ; 

P{ r = n } = (1 - P ) 2 F{S n = 0} + o(Q n ). 

Proof. Let m be such that n = 2m in the case of even n and n = 2m + 1 
otherwise. We have 



[r = n }= Y ^{S m = x,r = n}. 
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Since S m and the sequence {S m+ k — S m ,k > 1} are independent, the xth 
summand in the latter sum is equal to the product 

^{S m = x,t > m} 

xF{S m+k - S m 7^ —x for all k = 1, . . . , n - m - 1, S n - S m = -x}. 

The second probability here is equal to 

F{S n - S m+k ^ for all k = 1, . . . , n - m - 1, S n - S m = -x} 
= f{Sk ^ for all k — 1, . . . , n — m — 1, S n - m = —x} 
= F{r > n - m, S n - m = -x}, 

where Sk := £i + • • • + £fc := £ n _fe+i, and r is the first return time to zero 

of the random walk Sk- The random walk has the same distribution as 

Sk, so 

P{ r = n} = ¥ i S ™ = x,t> m}F{S n _ m = -x,t > n - m} 

x£Z d \0 

= ^2 ^{Sm = x,r > m}F{S n - m = -x, t > n-m} 

x&Z d \0 

= Si + S 2 , (11) 

where Si is the sum over \\x\\ < logn =: r n and S 2 is the sum over \\x\\ > 
logn. By the condition (jlj) and regular variation of Q n , 

Si < J2 nS m = x}F{S n „ m = -x} 

\\x\\ <log n 

< Q m Qn- m {2 log n) d = o(Q n ) as n -)> oo. (12) 
By Lemma El as n — > oo, 

S 2 = (1-p) W{Sm = x,T>m}lF{S n ^ m = -x} + o(Q n )} 

\\x\\>\ogn 

= {1-p) Y W{Sm = x,T>m}F{S n ^ m = -x} + o{Q n ). 

\\x\\>\ogn 

Repeating these arguments to the first multiple we obtain that 

S 2 = (l-p) 2 Y nS m = x}F{S n . m = -x} + o(Q n ). 

\\x\\ >log n 
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Taking also into account (1121) we finally obtain that 

£ 2 = (l-p) 2 P{S n = 0} + o(Q„). (13) 

Substituting (fT2"|) and (fT3"|) into (TTT|) . we arrive at the desired conclusion. 

The proof of Theorem H] is now immediate, since by Corollary [2] the se- 
quence Q n = k/C n satisfies all conditions of Lemma [61 for suitable k. The 
proof is complete. 

Remark 8. The same question may be addressed in the more general setting 
of matrix norming with the help of results by Griffin [8]. The key point is 
that the norming sequence \B n \ in his Theorem 6.4 is automatically regularly 
varying: we owe this comment to Phil Griffin, in a private communication. 
It is then easy to see that our result extends to this situation whenever d > 3 
or d = 2 and transience is assumed. 
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